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ABSTRACT 

We evaluate the extent to which newly detected cxoplanctary systems containing 
at least four planets adhere to a generalized Titius-Bode (TB) relation. We find 
that the majority of exoplanet systems in our sample adhere to the TB relation 
to a greater extent than the Solar System does, particularly those detected by the 

-—. Kepler mission. We use a generalized TB relation to make a list of predictions for the 

existence of 133 additional exoplanets in 65 multiple-exoplanet systems: 68 candidates 
from interpolation, 65 candidates from extrapolation. We predict the existence of 2 
low-radius (R < 2.3i?®) exoplanets in the habitable zone of KOI-490 and that the 
average number of planets in the habitable zone of a star is 1-2. The usefulness of 
the TB relation and its validation as a tool for predicting planets will be partially 

^h tested by upcoming Kepler data releases. 

o 

r~^ 1. INTRODUCTION 

During the last few years the number of multiple-exoplanet systems has increased rapidly. 

rG While most multiple-exoplanet systems contain two or three planets, a rapidly growing number 

II of systems now contain at least four planets. The architecture of these systems can begin to 

be compared to our Solar System. Many exoplanet systems appear to be much more compact 

■*— > and their planets more evenly spaced than the planets of our Solar System. This paper sets 

out to quantify these tendencies and place constraints on the location of undetected planets in 



> 



multiple-exoplanet systems. 



The approximately even logarithmic spacing between the planets of our Solar System mo- 
tivated the Titius-Bode (TB) relation, which played an important role in the discovery of the 
Asteroid Belt and Neptune, although Neptune was not as accurately represented by the TB re- 

i f J lation as the other planets (????). Uranus could have been discovered earlier if the TB relation 

had been take more seriously (?). Since the TB relation was a useful guide in the prediction 
of undetected planets in our Solar System, it may be useful for making predictions about the 

C^) periods and positions of exoplanets in multiple-exoplanet systems. The main goal of this paper is 

to test this idea. We first examine the degree to which multiple-exoplanet systems adhere to the 

£> TB relation. Due to limited detection sensitivity, we have no reason to believe that all exoplanets 

in a given system have been detected. After assessing this incompleteness in several ways, we 

find that the more complete the exoplanet system, the more it adheres to the TB relation. If this 
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trend is correct, the TB relation can be profitably applied to multi-planet systems to predict the 
periods, positions and upper mass or radius limits for as yet undetected exoplanets. 

The distribution and separation of planets in our Solar System can be parameterized by a 
generalized (two parameter) TB relation (??): 

a n = aC n , n = 0,1,2,3,... (1) 

where a n is the semi- major axis of the nth planet or satellite, a is the best-fit parameter associated 
with the semi-major axis ao of the first planet, and the powers of C parametrize the logarithmic 
spacing. This form of the TB relation involves two free parameters (a and C), instead of three 



from the original TB relation (?, Equation Al), which was derived emperically based on our 
Solar System only. 

The paper is organized as follows. Section [2] discusses the physics behind the TB relation. In 
Section [3] we review recent work on the TB relation and planet predictions that have been made 
based on it. In Section H] we describe our analysis and present our main results (predictions for 
cxoplanet periods). We discuss the details of several systems, caveats and compare our results 
with previous studies in Section [5j We summarize our results in Section [6} 



2. PHYSICS BEHIND THE TB RELATION 

As proto-planetary oligarchs accrete mass through collisions within a protoplanetary disk, 
they clear out the material in nearby orbits (as in the IAU definition of a planctFl), thereby 
excluding the presence of nearby planets. Starting from a relatively smooth disk (e.g. a sur- 
face density £ oc r~ 3 ' 2 ), the areas of mutual exclusion grow with time as the system relaxes 
dynamically. This produces a non-random distribution of planets in which the planet periods 
and positions are roughly logarithmically spaced (at least in the case of our Solar System, ?). 

A roughly logarithmic spacing between planet periods can be roughly parametrized by the 
generalized TB relation. Due to the stochastic nature of planetary scattering, planet migration, 
high eccentricities, and other variables of planetary formation, it is not obvious that the TB 
relation should even approximately fit exoplanet systems. 

From Equation 111 and Kepler's 3rd law (P n oc a n ), the generalized TB relation can be 
restated in terms of the periods of the orbits, 

P n =Pa n , n = 0,l,2,...N-l (2) 

where P n is the period of the nth planet or satellite, and P and a(— C 3 / 2 ) are parameters to be 
fit for each system. 

It is well known that the period ratios of planets and satellites in the Solar System show 
a preference towards near mean motion resonance (NMMRs) (???). More recently it has been 
shown that planets in exosystems also exhibit the same preference towards near mean motion 

1 http://www.iau.org/static/rcsolutions/Rcsolution_GA26- 5-6.pdf 



resonances (???). From Equation [2] the period ratio between adjacent orbits is P n+ i/P n = a. If 
a is approximately constant for all n (all planet pairs in the system), then the generalized TB 
relation will fit the system well. The physical mechanisms underlying these empirical observations 
are not fully understood. 

For example, Jupiter and Saturn are in a 5 : 2 NMMR. Approximately the same NMMR 
exists between Mercury and Venus, Mars and the Asteroid Belt and the Asteroid Belt and 
Jupiter. The remaining planet pairs have period ratios of ~ 1.6 (~ 3 : 2), ~ 1.9 (~ 2 : 1), ~ 2.9 
(~ 3 : 1) and ~ 2.0 (~ 2 : 1) for Venus and Earth, Earth and Mars, Saturn and Uranus and 
Uranus and Neptune respectively. The larger the scatter around the dominant NMMR, the less 
well the system adheres to the generalized TB relation. In the case of our Solar System, the 
smaller the scatter of period ratios around the dominant 5 : 2 ratio, the better the fit to the 
generalized TB relation. 



3. PREVIOUS WORK USING THE TB RELATION 

? analyzed adjacent planet period ratios in 11 simulated systems. Each system was evolved 
with a 1 Mq host star but with different planet mass distributions. He suggested that the TB 
relation results from some NMMRs (5 : 2 and 3 : 2) being more stable than others. Instability 
tends to exclude certain period ratios between interacting planets. The stronger the interactions, 
the more exclusion in period ratio space, the narrower the bunching of periods ratios, and 
the better the fit to the TB relation. ? performed simplified numerical simulations of 10 5 
planetesimals distributed in systems with different initial planetesimal surface densities S(a), 
where a is the semi-major axis. A merger resulted when two planetesimals underwent a close 
encounter and simulations continued until no more close encounters were possible. The final 
location of planets could be well fit by a TB relation when the surface density S(a) ex a~ 3//2 . 
This same form of S is derived for the Solar System from the Minimum Mass Solar Nebula (??) 
and a similar form, S(a) oc a~ 16 has been recently derived for a Minimum Mass Extrasolar 
Nebula from Kepler data (?). 

? created model planetary systems by injecting accretion nuclei into a disk of gas and dust 
until the dust was depleted. They found that simulated systems adhered to the generalized TB 
relation "about the same as the Solar System." Another numerical approach was taken by ? 
where model planetary systems were generated and subjected to simple stability criteria. For 
example, systems whose adjacent planets were separated by less than k times the sum of their 
Hill radii, where k ranged from to 8, were discarded. They found that the adherence to the TB 
relation increases with k, i.e. as the region of semi-major axis space open for adjacent planets is 
restricted, the adherence to the TB relation increases. In period space, as the effect of mutual 
exclusion increases, the NMMRs closer to 1 are excluded and this results in a narrowing of the 
spread of period ratios in the system. 

? looked for the possibility of additional planets in the five-planet 55 Cnc system based on 
the generalized TB relation (Equation [I]). They predicted planets at ~ 2 AU and at ~ 15 AU but 
made no attempt to quantify the uncertainties of those predictions. ? applied the generalized TB 
relation to multiple-exoplanet systems observed with HARPS and found reasonable fits, although 



they made no planet predictions. ? applied the three-parameter TB relation (Equation Al I to 
55 Cnc. He proposed that should the TB relation be applicable to the sytem, undetected planets 
may exist at semi-major axes of 0.085 AU, 0.41 AU, 1.50 AU and 2.95 AU. Additionally, a planet 
exterior to the outermost detected planet was predicted to exist between 10.9 and 12.6 AU. For 
a comparison of the two and three parameter TB relations see Section |5.5| and Appendix [XI 

? analyzed the period ratios of pairs of planets in individual multiple-planet systems, com- 
paring the distributions of the Solar System and 55 Cnc to that expected from a generalized 
TB relation (Equation [I]) . ?, repeated the analysis for 31 multiple-exoplanet systems. Even 
without taking into account any correction for the expected incompleteness of planet detections, 
? tentatively concluded that the adherence of the planets of our Solar System to a TB relation 
is not "due to chance" and that one cannot rule out the possibility that the TB relation can be 
applied to exosystems. 

By utilizing numerical N-body simulations one can estimate the regions of semi-major axis 
and eccentricity space which could host an additional planet without the system becoming dy- 
namically unstable. Although planets do not necessarily have to occupy these stable regions, ? 
argue that most planetary systems lie near instability, leading to the 'packed planetary systems' 
(PPS) hypothesis (????). The PPS hypothesis suggests that the majority of systems will be 
near the edge of dynamical instability. That is, undetected planets will occupy a stable region 
and will tend to bring the system closer to (but not over the edge of) instability. 

The PPS hypothesis and the TB relation are not orthogonal conceptions. As the spacing 
between planets decreases, i.e. as a system becomes more tightly packed, a system's adherence to 



the TB relation increases (see Section 4.2 1. The PPS hypothesis suggests that planets are spaced 
as tightly as possible without the system becoming dynamically unstable. Since dynamical 
stability is correlated with planets being evenly spaced, we may then expect systems which 
adhere to the PPS hypothesis to be more likely to adhere to the TB relation. However, the 
degree to which PPS systems adhere to the TB relation is poorly understood since PPS analyses 
predominately focus on two-planet systems (??). 

While numerical simulations can be useful for predicting stable regions, computing con- 
straints typically confine direct orbit integrations to two-planet systems. It becomes inefficient 
to integrate a large number of systems, each containing a large number of planets. There are 
now 65 systems containing four or more planets, and this number is growing rapidly due to the 
Kepler mission. The remainder of this paper discusses our tool for analyzing a large number 
of systems and constraining regions where we expect undetected planets to exist based on our 
generalized TB relation. 

Here we improve on previous TB relation investigations first by more carefully normalizing 
the % 2 /dof fit of the TB relation to our Solar System and then by fitting the TB relation to 
65 exoplanet systems which contain at least four planets (Figure pi . We choose this minimum 
number of planets as a compromise between, on the one hand, having enough systems to analyze 
without suffering too badly from small number statistics, and on the other hand, having enough 
planets in each system to fit 2 parameters and minimize incompleteness of the systems. Seven of 
the 65 systems were detected by the radial velocity method; Gl 876, mu Ara, ups And, 55 Cnc, 
Gl 581, HD 40307 and HD 10180. One system, HR 8799, was detected by direct imaging. The 



remaining 57 are candidate systems detected by the Kepler mission. Four of the Kepler systems 
have been confirmed by radial velocity observations (Kepler-11, Kepler-20, Kepler-33 and Kepler- 
62). The remaining 53 Kepler detected systems consist of planet candidates. However, ? conclude 
that < 2% of Kepler candidates in multiple planet systems are likely to be false positives. 



4. METHOD AND RESULTS 

Since the jostling during planetary system formation is a messy process that can include 
planetary scattering and migration, we expect varying degrees of adherence to the TB relation. 
Previous TB relation studies of exoplanetary systems (??) chose where to propose additional 
planets based on the solution which minimizes % 2 /dof. If our Solar System is treated like an ex- 
osystem, where limited detection sensitivity precludes the detection of some planets, this method 
does not necessarily recover the complete Solar System. We use the Solar System as a guide to 
how well we expect exosystems to adhere to the TB relation. In log space, the period form of 
the generalized TB relation (Equation [2]) is given by 

logP„ = logP + nloga = A + Bn, n = 0,1,2, ...,N- 1 (3) 

where A and B are parameters being fit to each system. The x 2 /dof value for a system containing 
N planets can then be calculated by 



X 2 (AB) _ 1 ^ ( (A + Bn) - log P, 



N-l , , a , -r, \ i n \ 2 

i (4) 

n=0 v n / 

where a n is the unknown uncertainty associated with the various reasons that the log of the 
period of the nth planet does not conform to the TB relation. The er„ are not the uncertainties 
in the precisions of measurements of logP„ because the deviations of logP„ from a good fit 
to the TB relation are dominated by the uncertainties in the complicated set of variables that 
determine the degree to which the NMMRs of adjacent planets bunch up around the same ratio. 
This is poorly understood (??). In our analysis, we choose uniform uncertainties (a n — a) and 
scale a for each system by the mean log period spacing between planets in each system, which 
we refer to as the sparseness or compactness of the system (see Equation [6] and Figure [4]). N — 2 
is the number of degrees of freedom (N planets - 2 fitted parameters (A and B)). 

Given the periods of the planets in an exoplanet system that contains 4 planets (e.g. 
Po, Pi, Pi-, P3), our strategy is to insert a new fifth planet (or more) between two adjacent planets 
and then fit the new system (e.g. Po, Pi, P2, P3, P4, ...) to the TB relation to see if the new x 2 /dof 
from Equation [4] is lower. The maximum mass or radius of the inserted planet is determined 



by the minimum signal to noise ratio of detected planets in the same system (Section 4.3). An 
example of this is shown in Figure [T] where we remove the Asteroid Belt and Uranus from the 
Solar System and then insert planets and compute the % 2 /dof fit to the TB relation. When the 
Asteroid Belt and Uranus are included in the fit to the Solar System, the % 2 /dof decreases from 
~ 3.7 to ~ 1.0, thus improving the fit to the TB relation. We do not increase the dof when we 
insert planets because the new planets are inserted on the best fit line and do not contribute to 
X 2 - However, inserting planets does increase the compactness and therefore, as shown in Figure 
4, does reduce the value of a. In some cases the insertion of an additional planet into a system 



results in a decrease of % 2 /dof which means that the new system with the inserted planet adheres 
to the TB relation better than the system without the insertion. 

Of the many possible ways to insert planets, we want to identify the way which improves 
the adherence to the TB relation the most. At the same time, we want to avoid inserting too 
many planets, since the lowest x 2 /dof does not necessarily correspond to the complete system 
(at least in the case of the Solar System). See for example, panels c and g of Figure [I] To 
quantify this, we define a variable 7, which is a measure of the fractional amount by which the 
X 2 /dof improves, divided by the number of inserted planets (i.e. the improvement in the x 2 /dof 
per inserted planet): 




(5) 

where xf an d Xf are the x 2 values from Equation VM before and after the insertion of rii ns planets. 
7 is calculated for each insertion solution and the solution with the highest 7 is chosen. Using 
7 recovers the actual Solar System when the input system does not include Uranus and the 
Asteroid Belt (Figure [l]). 



4.1. Isolating The Most Complete Systems 

To test the adherence of exoplanetary systems to the TB relation, we want to minimize the 
effects of incompleteness on the test. Thus we identify a more complete sample using a dynamic 
spacing criterion. After identifying a sample of exoplanet systems or subsets of exosystems 
which are likely to be more complete, we test their adherence to the TB relation. A simple 
approximation of the stability of a pair of planets is the dynamical spacing A (??), which is 
a measure of the separation of two planets in units of their mutual Hill radius (Equation Bl, 
Appendix IB]) . The dynamic spacing A is a function of the planet masses and hence dependent 
on our radius-mass conversion for transiting planets (R w 1.35M 038 , derived from 26 R < 8i?® 
exoplanets with mass and radius measurements). 

For small eccentricities and small mutual inclinations, ? found the stability of simulated 
systems increased exponentially with A. A disrupting close-encounter between planets was likely 
within 10 8 years for A < 10. The same A < 10 stability threshold was observed by ?. The A 
values for all adjacent pairs of exoplanets in our sample is shown in Figure |4~Tj There is a drop-off 
below A ss 10, in agreement with the simulations. Also, NMMR pairs dominate the A < 10 
pairs. Planet pairs with small dynamic spacing, A < 10, are less likely to be stable, but if they 
are stable, they are more likely to be NMMR pairs. We use this criterion to determine which 
pairs are more likely to be in complete systems where the insertion of an additional planet (with 
an assumed mass of 1 M ffi ) between the pair would result in a lower probability of stability: 
A < 10. Note that this assumed mass for the inserted planet plays a minor role, since the mass 
of the detected planet is always larger and dominates in the calculation of A (Equation Bl, 
Appendix [Bl) . 
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Fig. 1. — Fitting the TB relation to the planets of our Solar System when Uranus and the Asteroid 
Belt are not included. This figure shows the ability of the TB relation to predict where planets could be 
inserted to reconstruct the incompletely detected planetary system. We use Equation Wl to compute % 2 
with dof =7 — 2 = 5. We normalize the Solar System's adherence to the TB relation by adjusting a to 
ensure that x 2 /dof = 1 when Uranus and the Asteroid belt are included in the analysis (dof =9 — 2 = 7). 
The removal of Uranus and the Asteroid Belt increases the x 2 /dof from 1 to 3.7. The green and blue 
lines show the best fit generalized TB relation (minimizing \ 2 /Aoi of Equation WJ) before and after the 
insertion of 1 (top row), 2 (second row), 3 (third row) and 4 (fourth row) planets. We use the maximum 
value of 7 (Eq. 5) to prioritize these fits. The parameter a has been scaled in each panel according to 
the sparseness/compactness of the system (Equation 6) as described in Figure 4. This figure of the Solar 
System has been constructed to be as comparable as possible to our method of predicting the periods 
and positions of undetected exoplanets in multi-planet exoplanetary systems. Figure [2] shows the same 
method applied to the 55 Cnc system. 
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Fig. 2. — Illustration of our method for the 55 Cnc system, which has an initial x 2 /dof of 1.54. One to 
four planets are inserted between detected planets and new a and xVdof values are calculated for each 
insertion solution. For each number of planets inserted from 1 (top panels) to 4 (bottom panels), the two 
lowest x 2 /dof solutions are displayed. The insertion solutions with the two highest 7 values (Equation 
5) are indicated by a white background. An approximate habitable zone is represented by the horizontal 
green bar, which represents the effective temperature range between Venus and Mars, using an albedo 
of 0.3. 
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Fig. 3. — Dynamical spacing A for adjacent planet pairs within exoplanet systems containing four or 
more planets, separated into NMMR and non-NMMR planet pairs. The threshold to be in NMMR is 
set arbitrarily at x < 2%, where x = \Nj/Ni — P n+ i / P n \ / (Nj / Ni) . N and Nj are positive integers with 
Ni < Nj < 7 and P n and P n +i are the periods of the planet pair. There is a drop-off of planet pairs at 
A < 10 as found in simulations (??). The values for the Solar System are shown for reference. 



If the insertion of an additional planet between each pair in a system results in A < 10 for 
all pairs, then that system is likely to be complete without any insertions. Eleven systems satisfy 
this criterion, they are KOI-116, KOI-571, KOI-623, KOI-671, KOI-720, KOI-730, KOI-1358, 
KOI-2038, KOI-2169, KOI-2722 and HR 8799. A further 16 systems have 3 or more adjacent 
planet pairs with A < 10 (after insertion) and we treat these subsets (without insertions) as 
their own complete system, giving a total of 27 systems in our most-complete sample. 

Of the 27 systems in our most-complete sample, 22 fit the TB relation better than the Solar 
System (% 2 /dof < 1, with a scaled with compactness as shown in Figure 4). Two systems fit 
approximately the same as the Solar System and 3 of the 27 systems fit worse than the Solar 
System (% 2 /dof > 1.5). If additional planets exist in these systems, they are likely to be in a 
near mean motion resonance with one or more of the detected planets. 

Considering only the systems which are most likely to be complete, 89% adhere to the TB 
relation to approximately the same extent or better than the Solar System. Thus planetary 
systems, when sufficiently well sampled, have a strong tendency to fit the TB relation. Taking 
this strong-tendency-to-fit as a common feature of planetary systems, we make predictions about 
the periods and positions of exoplanets that have not yet been detected. Specifically, we fit the 
TB relation to systems that are less complete, insert planets into a variety of positions and find 
the positions that maximize the 7 value. Figure [T] illustrates our procedure for the Solar System 
when Uranus and the Asteroid Belt are missing. Figure [2] and Table [T] illustrate our procedure 
applied to the 55 Cnc system. The usefulness of planet predictions based on the TB relation 
depends on how much of a better fit each insertion solution produces (Eq. 5). 

To test the validity of this idea further, we removed individual planets from systems in our 
most-complete sample, thus creating less-complete systems. We then inserted a planet in the 
location predicted by our procedure of maximizing 7. The predicted location was between the 
correct pair of adjacent planets 100% of the time. The predicted period (a new P n in Equation]!]) 
was within our la of the original period ~ 79% of the time. 



4.2. The Compactness of Planetary Systems 

One factor not considered by other studies is the relation between the sparseness or com- 
pactness of a system and its adherence to the TB relation. Our initial analyzes of exoplanet 
systems used the value of a that made x 2 /dof = 1 for the Solar System. However, this led 
to x 2 /dof << 1, particularly for the most complete (and compact) exoplanetary systems. The 
Solar System is much more sparse than these systems (Figure 5). We have included the Solar 
System in Figure 5 to show this difference. In more compact systems with shorter periods, one 
expects the tendency towards commensurability of neighboring orbits to propagate more easily 
to other orbits (??). This is especially important when analyzing systems discovered by the 
Kepler mission, where individual systems have multiple planets within ~ 50 day periods (~ 0.3 
AU). Thus, we introduce a scaling for a dependent on how compact or sparse the system is. We 
define the sparseness S p of a system as, 



log Pjy-i - log P 

N 



S p = ICr ( 6 ) 



10 



where Pjv-i and Pq are the largest and smallest planet periods in the system respectively and 
N is the number of planets in the system. 

In Figure HI we use the Solar System as a standard and make a linear fit between it and 
the (0,0) point. For each system in our analysis, the a in Equation H is calculated from this 
linear fit (<r = 0.270^), given the sparseness/compactness of the system (Equation 6). More 
compact systems are assigned smaller a values. Systems above the line adhere to the TB relation 
less well than the Solar System does, while those below the line adhere to the TB relation more 
closely than the Solar System does. Each time a hypothetical planet is inserted between detected 
planets, the system becomes more compact and a new, lower a is calculated. The periods of 
the inserted planets are calculated from the best fit TB relation. This ensures that the inserted 
planets make no contribution to the x 2 value. 

If a system has a % 2 /dof similar to or lower than the Solar System (< 1) before any insertions, 
no insertions arc made. Applying our 7 method to our sample results in insertions being made 
in 26 out of 65 exosystems. For all systems in our sample we predict the location of the next 
outermost planet (extrapolation). Our planet predictions are shown in Figures p] & \& and Tables 
2& 3. 



4.3. Upper Mass or Radius Limit 

We can place constraints on the maximum mass (radial velocity detections) or maximum 
radius (transit detections) of our predicted planets. M max and R max are calculated by applying 
the lowest signal-to-noise ratio (SNR) of the detected planets in the same system, to the period 
of the inserted planet. That is, we calculate the maximum mass or radius of a planet at the 
predicted period that could avoid detection based on the lowest SNR of the detected planets in 
that same system. For transiting planets in the same system, SNR ex r 2 P -1 / 2 where r is the 
planet radius and P is the planet period. For radial velocity detected planets in the same system, 
SNR oc KP~ 5 ' e oc mP~ 7 ' e where K is the semi-amplitude of the radial velocity and m is the 
planet mass. Maximum masses and radii for predicted planets can be found in Table 2. 
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0.4 
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Fig. 4. — Top Panel: a used in Equation Rl as a function of the average log period spacing between 
planets, S p (Equation Rjl, of multiple planet systems. The thick line goes through two points: one point 
is the origin (0, 0). The other point (S p , a), is the sparseness of our Solar System (Equation™ and the 
a required for our Solar System to yield xVdof = 1 in Equation Ul The Solar System is shown as a blue 
triangle while the four small blue squares above it (from top to bottom) are the values of (S p ,a) from 
the Solar System if Saturn, Jupiter, Uranus and the Asteroid Belt are individually removed. The main 
prograde satellites of Jupiter, Saturn, Uranus and Neptune are also fit and indicated by small green 
crosses and labeled " Jup" , "Sat" , "Ura" and "Nep" respectively. The systems in our most complete 



sample (see Section 4.1l are indicated by green dots. 



Bottom Panel: Similar to the top panel but with the following exceptions. All systems are shown rather 
than only our most complete set. Black dots indicate systems where no insertion was made. Gray dots 
indicate a system where an insertion was made, and red dots indicate the (S p ,a) of the system after 
insertions have been made. Each gray point is connected to a red point by a gray arrow. The dashed 
box represents the range of the top panel, which encompasses all points in our most complete sample. 
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Fig. 5. — Orbital periods of planets in multiple planet systems containing at least 4 detected planets. 
Systems where planet insertions have been made are sorted in descending order of the highest 7 found in 
each system (Equation [5] and Table 2). Systems without planet insertions (lower 2/3 of figure) are sorted 
in ascending order of x 2 /dof (Equation 4). Inserted and extrapolated planets are shown with red filled 
circles and red open circles respectively. The 7 value for the Solar System is calculated by excluding, 
then including, the Asteroid Belt. "*" indicates at least two adjacent planet pairs in the system have 
A values < 10 if we had inserted a planet between each pair. "**" indicates all adjacent planet pairs 
have A values < 10 if we had inserted an additional planet between each pair (see Section E|. For the 
systems where no planet was inserted by interpolation, predicted extrapolated planets are listed in Table 
3. Planet letters are shown for Gliese 876. 

1 The two middle planets in KOF2248 (bottom of plot) are dynamically unstable. One or both of the 
candidates may not be a planet (?). 
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Planet T Eff (K) 



Fig. 6. — Same as the previous figure except that the periods of the planets have been converted 
into effective planet temperatures (assuming an Earth-like albedo of 0.3), based on the luminosity of 
the stellar host. The vertical green bar represents an approximate habitable zone using the effective 
temperature range between Venus and Mars, using an albedo of 0.3. We estimate that the average 
number of planets in the habitable zone (HZ) of a star is approximately 1-2. This estimate is based on 
i) the assumption that planetary systems extend across the HZ, and ii) on the number of HZ planets in 
the planetary systems shown here that span the HZ. 
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Tabic 2: Systems With Interpolated and Extrapolated Planet Predictions 



System 


Number 


7 


A-y a 


/x 2 \ 


/x 2 \ 


Inserted 


Period 


A 


R b 


M majc b 


T e ff 


Inserted 


(Eq. 5) 


U°*7i 


\d t) f 


Planet # 


(days) 


(AU) 


(a®) 


( M ffi) 


(K) 


KOI-1052 


2 


234.0 


1528.2 


1.36 


0.01 


1 
2 




11 ± 2 
28 ± 3 


0.10 
0.18 


1.6 

2.0 


- 


909 

669 














3 


E c 


108 ± 12 


0.46 


2.8 


- 


423 


Glicsc 876 


3 


210.7 


32.5 


7.88 


0.02 


1 
2 
3 




3.9 ± 0.7 

8 ± 2 
16 ± 3 


0.03 
0.05 
0.08 


- 


0.7 
1.2 
2.2 


489 
:iss 
308 














4 


E 


245 ± 40 


0.52 


- 


21.9 


122 


KOI-1952 


2 


171.1 


10.9 


3.26 


0.01 


1 
2 




13 ± 2 
19 ± 2 


0.10 
0.14 


1.5 
1.6 


- 


828 
720 














3 


E 


65 ± 7 


0.31 


2.2 


- 


474 


Kcplcr-62 


6 


107.3 


40.6 


4.04 


0.01 


1 
2 
3 
4 
5 
6 




8.5 ± 0.8 
27 ± 3 
40 ±4 
58 ± 6 
85 ± 8 

182 ± 18 


0.07 
0.15 
0.20 
0.26 
0.33 
0.55 


0.5 
0.6 

0.7 
0.8 
0.9 
1.0 


- 


649 
442 
389 
342 
301 
233 














7 E 


391 ± 37 


0.92 


1.2 


- 


180 


KOI-490 


5 


79.1 


21.7 


7.56 


0.02 


1 

2 
3 

4 

5 




13 ± 2 

39 ± 5 

66 ± 9 

114 ± 15 

197 ± 26 


0.10 
0.20 
0.29 
0.41 
0.59 


1.2 
1.5 
1.7 

2.0 
2.3 


- 


585 
406 
339 
282 
235 














6 


E 


584 ± 77 


1.22 


3.0 


- 


163 


KOI-248 


1 


18.3 


17.5 


2.48 


0.13 


1 




4.3 ± 0.5 


0.04 


1.4 


- 


634 














2 


E 


31 ±4 


0.16 


2.2 


- 


330 


KOI-500 


2 


15.2 


4.6 


5.42 


0.18 


1 
2 




1.5 ± 0.2 
2.2 ± 0.2 


0.02 
0.03 


1.1 
1.2 


- 


1055 

929 














3 


E 


15 ± 2 


0.11 


2.0 


- 


490 


KOI-1567 


2 


11.8 


12.1 


1.51 


0.07 


1 
2 




12 ± 2 
29 ± 3 


0.10 
0.18 


2.0 
2.5 


- 


66s 
494 














3 


E 


70 ± 8 


0.32 


3.1 


- 


366 


KOI-1198 


6 


10.8 


1.5 


7.19 


0.11 


1 
2 
3 
4 
5 
6 




1.5 ± 0.2 
2.3 ± 0.3 
3.3 ± 0.4 
4.9 ± 0.5 
7.3 ± 0.7 
24 ± 3 


0.03 
0.04 
0.05 
0.06 
0.08 
0.17 


1.2 
1.3 
1.4 
1.6 

1.7 
2.3 


- 


1855 
1626 
1425 
1249 
1095 
737 














7 E 


53 ± 6 


0.29 


2.8 


- 


566 


KOI-1306 


2 


8.8 


0.6 


4.14 


0.23 


1 

2 




10 ± 2 
18 ± 3 


0.09 
0.13 


1.4 
1.6 


- 


841 
700 














3 


E 


52 ± 7 


0.27 


2.1 


- 


485 


ups And 


2 


8.5 


1.0 


5.31 


0.30 


1 

2 




18 ± 6 

70 ± 22 


0.14 
0.36 


- 


3.8 

11.9 


847 
537 














3 


E 


16300 ± 5000 


13.57 


- 


1116.0 


87 


Kepler-20 


1 


7.4 


1.0 


3.51 


0.42 


1 




40 ± 6 


0.22 


1.2 


- 


500 














2 


E 


133 ± 19 


0.49 


1.6 


- 


332 


mu Ara 


4 


6.0 


1.8 


4.15 


0.17 


1 
2 
3 

4 




23 ± 5 

54 ± 11 

127 ± 26 

1690 ± 350 


0.16 
0.29 
0.52 
2.90 


- 


7.7 
15.8 
32.4 

280.2 


705 
529 

396 
167 














5 


E 


9500 ± 2000 


9.17 


- 


1180.7 


94 


Gl 581 


1 


5.0 


0.9 


3.72 


0.63 


i 




30 ± 5 


0.13 


- 


3.0 


238 














2 


E 


139 ± 23 


0.35 


- 


11.2 


141 


KOI-505 


3 


4.7 


0.6 


8.20 


0.56 


1 
2 
3 




22 ± 3 
35 ±4 
56 ± 7 


0.15 
0.20 
0.27 


4.6 
5.1 
5.8 


- 


788 
676 
580 














4 


E 


139 ± 16 


0.51 


7.3 


- 


426 


KOI-880 


1 


4.0 


2.4 


1.35 


0.28 


1 




12 ± 2 


0.10 


2.1 


- 


761 














2 


E 


117 ± 20 


0.45 


3.7 


- 


354 


Kcplcr-11 


1 


3.7 


1.7 


3.15 


0.69 


1 




75 ± 8 


0.34 


2.9 


- 


449 














2 


E 


171 ± 17 


0.59 


3.5 


- 


341 


KOI-1831 


3 


3.5 


1.7 


2.54 


0.23 


1 
2 
3 




6.5 ± 0.7 

9.8 ± 1.0 

23 ± 3 


0.06 
0.08 
0.15 


0.9 
1.0 

1.2 


- 


800 
697 
529 














4 


E 


78 ± 8 


0.34 


1.6 


- 


350 


Kcplcr-33 


1 


3.2 


0.1 


3.69 


0.90 


1 




8.9 ± 0.8 


0.09 


1.9 


- 


1162 














2 


E 


68 ± 7 


0.35 


3.2 


- 


590 


KOI-250 


5 


3.1 


2.5 


2.26 


0.14 


1 
2 
3 

4 

5 




4.9 ± 0.4 
6.7 ± 0.6 

9.2 ± 0.8 
25 ± 2 
34 ± 3 


0.05 
0.06 
0.07 
0.14 
0.17 


1.2 
1.3 
1.4 
1.8 
2.0 


" 


687 

617 
554 
401 
360 














6 


E 


64 ± 5 


0.27 


2.3 


- 


290 


KOI-245 


3 


2.9 


2.4 


1.56 


0.17 


1 
2 
3 




16.8 ± 0.9 
27 ± 2 
33 ± 2 


0.12 
0.16 
0.18 


0.3 
0.4 
0.4 


- 


582 
503 
467 














4 


E 


64 ±4 


0.28 


0.5 


- 


374 


55 Cnc 


2 


2.6 


0.5 


1.49 


0.25 


1 
2 




4 ± 2 
1080 ± 370 


0.04 
1.98 


- 


1.4 
178.2 


1117 

158 














3 


E 


20100 ± 6900 


13.97 


- 


2046.5 


59 


KOI-1336 


2 


2.5 


11.0 


1.07 


0.19 


1 
2 




6.8 ± 0.7 
26 ± 3 


0.07 
0.17 


1.7 
2.4 


- 


1053 

679 














3 


E 


61 ± 6 


0.31 


3.0 


- 


507 


KOI-952 


1 


2.2 


0.1 


2.36 


0.76 


1 




1.5 ± 0.3 


0.02 


0.9 


- 


902 














2 


E 


41 ± 7 


0.19 


2.1 


- 


298 


HD 40307 


4 


1.3 


0.2 


1.91 


0.32 


1 
2 
3 

4 




6.3 ± 0.7 
15 ± 2 
81 ± 9 

123 ± 13 


0.06 
0.11 

0.33 
0.44 


- 


0.4 
0.8 
3.4 

4.8 


814 
613 
348 
302 














5 


E 


287 ± 30 


0.78 


- 


9.7 


227 


KOI-719 


3 


1.2 


0.1 


1.32 


0.31 15 


1 
2 
3 




6.2 ± 0.6 
14 ± 2 
20 ± 2 


0.06 
0.10 
0.13 


0.6 

0.7 
0.8 


- 


693 
532 

467 














4 


E 


66 ± 7 


0.28 


1.1 


- 


315 





a A 7 = (71 — 72J/72 where "fi and 72 arc the highest and second highest 7 values for that system respectively. 

-Rmax and M max are calculated by applying the lowest SNR of the detected planets in the system to the period of the inserted planet. 

i.e. SNR oc r P~ ' for the planetary radius, where r is the planet radius and P is the planet period. SNR. oc KP~ ' for planets 

detected by radial velocity, where K is the semi-amplitude velocity (sec Section 4.3). 
c A planet number followed by "E" indicates the planet is extrapolated (has a larger period than the outermost detected planet in the 

system). 



Table 3: Systems With Only Extrapolated Planet Predictions 

Svstem (*") Period A R — a Mml " a Teff 

y \drt) t (days) (AU) (Re) (Me) (K) 

522 
374 
224 
413 
397 
349 
494 
293 
477 
468 
437 
640 

3730.5 35 

602 
425 
592 
508 
692 
461 
455 
336 
440 
499 
541 
613 
491 
381 

147.9 106 

630 
351 
868 
388 
309 
252 
408 
391 
865 
301 
565 

fi m „ and M max are calculated by applying the lowest SNR of the detected planets in 
the system to the period of the inserted planet, i.e. SNR <x r 2 p-!/ 2 for the planetary 
radius, where r is the planet radius and P is the planet period. SNR <x KP~ 5 ^ 6 
for planets detected by radial velocity, where K is the semi-amplitude velocity (see 
Section 4.3). 



5. DISCUSSION 

As exoplanet detection sensitivities improve, the number of multiple-planet systems will 
increase, along with the average number of planets detected in those system. Finding additional 
planets in these systems can be aided by making constrained predictions about their locations. 
This helps both in the observation stage where the sampling can be optimized to a specific period 
and in the analysis stage where more care can be taken at predicted locations. 

Some planet predictions have already been made by numerical integration. Regions are 
mapped in semi-major axis and eccentricity space where additional planets can be inserted while 
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0.11 
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0.12 


42 ±4 


0.19 
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0.14 


35 ±4 
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0.19 
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0.55 
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0.22 
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0.20 
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0.24 
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- 
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0.32 
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0.41 
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0.20 
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0.17 
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0.51 
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0.29 
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0.52 
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0.35 


2.0 
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228 ± 37 


0.77 
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0.56 


83 ±9 


0.37 


2.4 


[-1557 


0.58 


19 ±2 


0.12 


2.0 


[-1930 


0.60 


72 ±7 


0.35 


2.3 


[-671 


0.62 


27 ±2 


0.17 


1.5 


[-1563 


0.62 


27 ±3 


0.17 


3.7 


[-152 


0.64 


160 ± 16 


0.60 


5.1 


10180 


0.65 


6200 ± 1500 


6.67 


- 


[-1102 


0.75 


33 ±3 


0.20 


2.8 
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0.84 


117±17 


0.47 


2.5 


[-2169 


0.87 


7.7 ±0.4 


0.07 
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0.91 
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0.91 


250 ± 41 
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0.14 
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the system remains stable. Due to computing constraints, numerical integrations tend to focus on 
two-planet systems. For example, in our sample of systems containing four or more planets, only 
two systems have been analyzed for the stability of additional planets by numerical integration. 
As the number of planets detected per system and the number of systems continues to increase, 
directly integrating each system becomes more impractical. We emphasize that we do not expect 
the generalized TB relation can be applied to every type system (e.g. where exact resonances 
are present). 



5.1. Comparison With TB Relation Studies 

55 Cnc: 
? represented the semi-major axes of the planets in 55 Cnc with an exponential function (Equa- 
tion 



A3 1. They claimed a good fit when the 5th planet in the system was assigned n — 6, leaving 
a gap at n = 5 for an additional planet beyond the detection threshold. The fit was extrapolated 
to n — 7 to predict a planet beyond the outermost detected planet. The semi-major axes of the 
predicted planets were ~ 2 AU (P=l,086 days) and ~ 15 AU (P=22,309 days), with no attempt 
made to estimate uncertainties. After their paper was published, the period of the innermost 
planet was updated from 2.8 days to 0.7 days (?). 

When using the old value for the period, we similarly predict a single planet located at 
2.0 ± 0.5 AU (P=l,lll±406 days). When using the new value for the period, we predict two 
planets, the first at 0.04 ± 0.01 AU (P=3.1±l.l days), while the prediction for the planet at 
2.0 ± 0.5 AU (P=l,lll±406 days) is unchanged. We calculate the location of a planet exterior 
to the outermost detected planet at 15.1 ± 3.5 AU (P=22,988±7,816 days) and 14.7 ± 3.4 AU 
(P=22,079±7,492 days) when the old and new data are used respectively. All predictions are 
compatible with those made by ?, with the exception of our additional planet prediction at 
0.04 ± 0.01 AU (P=3.1±l.l days) caused by the updated period of the innermost planet. 



? applied the 3 parameter TB relation (Equation Al ) to 55 Cnc. The residuals were 
minimized when A ~ 0.037, B ~ 0.048 and Z ~ 1.98 and for planet numbers n = — oo, 1, 2, 4 
and 7. This implies undetected planets occupying the n — 0, 3, 5 and 6 spaces, calculated to 
have semi-major axes of 0.085 (P=9.5 days), 0.41 (P=100 days), 1.50 (P=705 days) and 2.95 AU 
(P=l,945 days). Cuntz predicted a planet exterior to the outermost planet (n — 8) located at 
11.8 ± 0.9 AU (P=15,599±l,780 days). The predicted planets at n = 3 and n = 8 are within the 
uncertainties of our predictions. The remaining three predicted planets are incompatible with 
our predictions, although no uncertainty estimate is made by Cuntz. 



5.2. Comparison With Numerical Simulations 

55 Cnc: 
Our prediction of a planet at 2.0 ± 0.5 AU (P=l,lll±406 days) is compatible with numerical 
simulations which have found that the 55 Cnc system can support an additional planet between 
0.9 (P=327 days) and 3.8 AU (P=2,844 days) with an eccentricity below 0.4 (??). ? also demon- 
strate that a planet exterior to the outermost detected planet will likely reside beyond 10 AU 
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(P=12,143 days), compatible with our more explicit prediction of 14.7±3.4 AU (P=22,079±7,492 
days). 

Gl 876: 
Numerical simulations by ? show a small stable region for an additional planet centered on 0.08 
AU (P=14.3 days) around the M dwarf Gliese 876. This region is stable due to the 2 : 1 mean 
motion resonance with Gliese 876 c, the <~ 0.7 Jupiter mass planet located at 0.13 AU (P=30 
days). ? initially found a similar stable region around 0.08 AU (P=14 days). Further analysis 
suggested that initially stable orbits became chaotic and unstable over larger time frames. In our 
analysis, we predict 3 additional planets between Gliese 876 c and the innermost planet, Gliese 
876 d. The inserted planets have semi-major axes ranging from 0.03 (P=3.2 days) to 0.1 AU 
(P=20 days). Although the numerical simulations of ? suggest that these planets are located 
within an unstable region, our predictions correlate to their most stable areas within that zone. 



5.3. Potential Habitable Zone Predictions 

KOI-490: We predict two planets located in the habitable zone of KOI-490. Ignoring 
systematic uncertainties, the non-detection of the proposed planets requires that they have radii 
< 2.0 i?^ and < 2.3 i?® with effective temperatures of ~ 280 K and ~ 235 K respectively. This 
makes KOI-490 an interesting system for the upcoming Kepler data releases. 

Extrapolated Systems: For systems which adhered to the TB relation well without any 
insertions, we predicted the period of the next outermost planet in the system. This resulted in 
a prediction of a potentially habitable terrestrial planet in three Kepler candidate systems. KOI- 
812 is a four-planet system with x 2 /dof = 0.04 (compared to 1.0 for the Solar System). When 
we insert a 5th planet into the system, it has an estimated effective temperature (assuming an 
Earth-like albedo of 0.3) of ~ 224 K and a maximum radius of 2.5 R© based on its current 
non-detection. KOI-571 (x 2 /dof = 0.12) is another four-planet Kepler candidate system where 
we have predicted an extrapolated planet. In this case, we predict a 5th planet with an effective 
temperature of ~ 293 K and a maximum radius of 1.8 R®, in the habitable zone of KOI-571. 
KOI-904 is a five-planet Kepler candidate system. We predict an extrapolated planet with an 
effective temperature of ~ 252 K and a maximum radius of 2.2 R©. 



5.4. Kepler Candidates With The Highest Probability Of Validation 

Assuming a high degree of coplanarity (?), and using the angular- momentum- weighted av- 
erage inclination of a planets in each system, <i>^, we compute a critical semi-major axis a cr it- 
Beyond this value, a coplanar planet would no longer transit (impact parameter b > 1, where 
b = acos<i>L/R* and i?* is the estimated radius of the host star). We compute these a cr it values 
for each transiting system in our sample. Kepler-62, KOI-1052, Kepler-33, Kepler-11 and KOI- 
1198 are the Kepler candidate systems where an undetected planet exterior to each outermost 
detected planet is most likely to transit i.e., the semi- major axis a out of the outermost detected 
planet contributes to a low value of a ut/flcrit- Similarly, extrapolated planets least likely to 
be seen by Kepler are in the following systems: KOI-2248, KOI-2169, KOI-904, KOI-952 and 
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KOI-248. 



5.5. The Two and Three Parameter TB Relations 

To check the robustness of our 2 parameter TB relation predictions, we repeated our analysis 
using the 3 parameter TB relation (Equation IaI]). The two forms of the TB relation do not yield 
identical results. For example, when the 3 parameter TB relation is applied to the Gliese 876 
system, no planet insertions are made. This is a result of the decoupling of the distance to the 
first planet from the rest of the system (Figure [71 Appendix [A} . The 2nd, 3rd and 4th planets 
in Gliese 876 fit the 3 parameter TB relation well while there is a large gap between the 1st 
and 2nd planets. However, this does not affect the fit since the distance to the first planet is 
its own decoupled parameter. This is not unique to the Gliese 876 system and the 3 parameter 
TB relation often does not insert a planet between the first and second detected planets, even if 
a large gap is present. Of the 25 systems where we predicted interpolated planets, 16 of these 
systems are also predicted to have interpolated planets by the 3 parameter TB relation. In 6 
of the 9 systems where planets are inserted by the 2 parameter TB relation but not by the 3 
parameter TB relation, there is a significant gap between the first and second planet in each 
system. The prioritization of the lists of predicted planets by the 2 and 3 parameter TB fits also 
varies due to the different 7 values produced for a given system by each form of the TB relation. 



6. CONCLUSION 

We first identified a sample of exoplanet systems most likely to be complete based on the 
dynamical spacing A of planet pairs. We applied the TB relation to each system in the sample 
and found that ~ 89% of the systems adhered to the TB relation to approximately the same 
extent or better than the Solar System. This was taken as evidence that the TB relation can be 
applied to exoplanet systems in general. 

We then applied the TB relation to the systems most likely to be incompletely sampled and 
we made planet predictions in these systems. Our method involved inserting up to 10 hypothetical 
planets into each system and we chose the insertion solution which maximized 7 (Equation pj , 
which is a measure of the improvement in the x 2 /dof value (Equation B per planet inserted. We 
made planet predictions (interpolations) inserting 68 planets into 26 out of 65 systems. For all 
systems we predicted the period of the next outermost planet (extrapolation) . 

The maximum mass or radius for each predicted planet was calculated based on the detection 
limit in the host system. Since the TB relation may be mass-independent (e.g. the Asteroid 
Belt), our predicted planets do not have a lower mass limit. Thus, some significant fraction of 
our predictions probably correspond to low mass planets or Asteroid Belt analogues, in which 
case, validation in the near future is unlikely. Our predictions include a number of potentially 
habitable terrestrial planets. We identified Kepler candidate systems where an undetected planet 
exterior to the outermost detected planet has the greatest chance of transiting, assuming near 
coplanarity. 
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As the number of candidate planets in a system increases, the system's light curve or radial 
velocity curve becomes more complicated. This increases the difficulty of detecting additional 
planets within the system. Using the TB relation to predict planetary orbital periods (Tables 2 & 
3) narrows the region of parameter space to be searched, and increases the likelihood of detecting 
an additional planet. We expect our predictions to be partially tested using new data from the 
Kepler mission in the near future. If a significant number of our predictions are substantiated, 
we expect the TB relation to become an important detection tool for systems which already 
contain a number of planets. 



A. THE TWO AND THREE PARAMETER TB RELATIONS 

The original TB relation was given a mathematical form by ?, 

a n = A + BZ n , n = -oo,0, 1,2,... (Al) 

where a n is the semi-major axis of the planet. The index n — — oo is assigned to the first planet, 
n = to the second planet, n = 1 to the third planet and so on. Historically, A — 0.4, B = 0.3 
and Z = 2.0 for the Solar System, while the best fit values are A ~ 0.382, B ~ 0.334, Z ~ 1.925. 
The two parameter TB relation can be written as 

a n = aC n , n = 1,2,3, ... (A2) 

The exponential function given by ? is equivalent to Eq. |A2| 

a n = ae Xn , = aC n n = 1, 2, 3, ... (C = e A ) (A3) 

By comparing Equation |A1| and Equation |A2| we can see they are most similar for large n. 



Therefore the extra term in Equation AT and the numbering methodology n — — oo,0, 1,... 
produces a difference in the fits to the inner part of a planetary system. This has led to the 
location of the Earth to be called peculiar (?). 

Other analyzes of the Solar System tend to use the generalized 2 parameter TB relation as 
we have done in this paper (?????). 

aC(C-l) aC 3 (C-l) ^^ A B(Z-l) BZ 2 (Z-1) 

-o < e o e — e e 



a(C-l) aC 2 (C-l) B BZ(Z-l) 



(a) The 2 parameter TB relation (Equation A2 i (b) The 3 parameter TB relation (Equation 



All 



Fig. 7. — Illustration of the parameters in the two forms of the Titius-Bode relation (Equations Al 
and |A2[ |. Labels show the distances between adjacent objects. The large filled circle represents the star 
while the smaller open circles represent the planets. 
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B. THE DYNAMICAL SPACING A 

The dynamical spacing A (??) between two planets in the same planetary system with 
semi-major axes a\ and ai is, 

A=^p (Bl) 

where Rh is the two planets' mutual Hill radius, given by 

where mi and 777,2 are the masses the inner and outer planets respectively. M* is the mass of the 
host star. 
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